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Abstract. The Jacobi theta-functions admit a definition through the autonomous diffe¬ 
rential equations (dynamical system); not only through the famous Fourier theta-series. 

We study this system in the framework of Hamiltonian dynamics and find corresponding 
Poisson brackets. Availability of these ingredients allows us to state the problem of a cano¬ 
nical quantization to these equations and disclose some important problems. In a particular 
case the problem is completely solvable in the sense that spectrum of the Hamiltonian can 
be found. The spectrum is continuous, has a band structure with infinite number of lacunae, 
and is determined by the Mathieu equation: the Schrodinger equation with a periodic cos- 
type potential. 
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trum of Hamiltonian 
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Why do we attach importance to a possi¬ 
bility to describe the classicallimit through 
theta-functions? .. . 

Thus, when quantizing integrable models, 
sooner or later we should come to a quan¬ 
tization of theta-functions [15, § 1]. 


1 Introduction 


Presently, the only known definition and introductory motivation to the famous theta-functions 
are through the corresponding Fourier series 


6iz\T) 


E 7ri(fc^T+2fc2) 


—oo 


Such series appear in numerous areas of mathematics and physics as building blocks for repre¬ 
sentation of solutions and also have a great interest in their own rights because of rich structural 
and applied properties. To be precise, there are some mathematical arguments to the series 
structure which come from considering the group and functional properties of Abelian varieties. 
The displayed Fourier structure then becomes, in a certain sense, unique and is generalized 
to multidimensional cases. In this work we shall deal only with Jacobian, i.e., 1-dimensional 
theta-functions. 

There is one imperfection of these ‘serial’ ways of introducing the 0-functions however. These 
definitions are pure functional, not invariant, and do not even look like a definition for any special 
function of mathematical physics. The latter is determined usually by a certain differential 
equation whereas the 0-series are not; e.g., they appear (always via theta-ratios) to represent 
solutions of many integrable models. On the other hand, we know that property of the model 

*This paper is a contribution to the Special Issue on Algebraic Methods in Dynamical Systems. The full 
collection is available at http://www.emis.de/journals/SIGMA/AMDS2014.html 
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to be algebraically integrable [17] is an invariant property (commuting operators) and is of 
a radically different kind from those models which are solvable in terms of various special 
functions [4], 

A further aspect is the fact that 0-functions are the nearest generalization of the elementary 
ones - harmonic sine/cosine-functions - which are related to the simplest classical solvable model 
and extensively used, through the harmonic quantum prototype 

= E'i', 

in the quantum field theories. Their low energy limits and excitations are treated as quantum 
mechanics with the field boson particles: statements of this elementary model. 

Yet another facet is a relation with 0-series quotients. It is known that basic 0-quotients are 
proportional to Jacobi’s elliptic sn-, cn-, and dn-functions [1, 11, 18] 

01 02 03 1 

— ~ sn, ~ cn, ~ ~ an, 

(74 (74 (74 

and describe dynamics of some famous models in the rigid body theories [9, 11]. In turn, most 
popular of them - Eulerian top - may be quantized, thereby determining the space quantization 
phenomenon, so called quantum rotator [10]. 

These and many other aspects give rise a problem towards the direct quantization of theta- 
functions themselves [15] if it is possible. Fortunately this is so indeed and we shall exhibit 
a particular case that is completely solvable in the quantum sense: spectrum of a quantum 
Hamiltonian can be found. The starting point here is the fact that 0-series posses the nice 
differential properties like elliptic sn-, cn-, dn-functions mentioned above. In the context above 
one can even say that any kind of differential definition to the 0-functions might be considered 
as more appropriate way of their introduction at all. 


2 Definitions and base dynamical systems 


Complete differential properties of Jacobian series were described very recently [5] and we re¬ 
produce them in a nutshell. These four 0-series are defined as [1, 18] 
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and satisfy the system of ordinary differential equations 
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where O'l is an additional (fifth) independent function and (?9, rf) are so called iJ-constants. In 
standard theory these ‘constants’ are the values of 0 (z|r)-functions at z = 0, i.e., i^k = ^fc( 0 |T), 
but for the differential formulation ( 1 ) such a meaning need not be implied; all the {r],9k) are 
allowed to be free external parameters [5]. In the abstract context the differential introduction 
of functions is much more preferable because we may consider these equations merely as a vector 
field and freely change coordinates on a manifold where the field acts. 

For a quantization procedure it is desirable to simplify vector fields into the polynomial ones. 
It appears that the sought-for transformation ( 01 , 02 )^ 3 ) 04,^0 ^ z,u,^) does exist and 

has the form 


|x = - 7 rd 3 i? 4 ^, y = - 7 r- 02 i? 4 ^, = -vrd 2 ' 03 ^, u = 0i, ^ + Atj 

under certain value of the parameter A(i?, 77 ). Let dot above symbol stand for a t-derivative. 
Then this substitution brings the 0-system above into a very simple polynomial form 

X = yz, ij = xz, z = xy, ^ u = ( 2 ) 

which has an embedded structure like ‘Russian dollies’ {{{x,y,z),£,),u). Moreover it is not 
difficult to show that the first three equations in ( 2 ) are nothing but the normalized (parameter- 
free) model of a free Euler’s rigid body [9, II] 



with the standard Hamilton and Casimir functions 

1/^2 Y‘^ Z‘^\ „ . „ 

H = - { — + — + — , K = 

2 \ A B C )' 

and the well-known Poisson bracket 
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The transformation (x, y, z) ^ (X, T, Z) has a linear form 

A^BC{C-B) B^CA{A-C) 

^{C-B){A-C){B-Af' ^{C-B){A-C){B-Af' 

^ _ Cy'ABjB - A) 

^{C-B){A-C){B-Af 

and holds only for the full asymmetric top parameters A ^ B ^ C ^ A. Quantization of 
this dynamical system is very well known as the quantum rotator [10]. We thus arrive at 
a quantization problem related to a 5D-extension of the classical 3D quantized rotator. 

As far as we know, no regular algorithm for searching for a Poisson bracket exists. Even 
Euler’s top is not an exception since its Lie algebraic SO(3)-bracket is ad hoc; from the Lie 
algebra 5 o( 3 ) to the Poisson structure (3). On the other hand, it is naturally to expect that 
knowing the solution in full details one could be able to find out all the associated constructions. 
Particularly, classical and all the other Poisson structures because they are closely related to 
integrals of dynamical systems. We thus need a complete set of these functions. 
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3 Integrals of motion 

Since equations (2) is an odd dimensional system we restrict ourselves in this work to its 4-di¬ 
mensional subsystem 

X = yz, ij = xz, z = xy, ^ = —x^. (4) 

The point here is that there are some serious foundations that even the full 5D system (2) may 
not be considered as closed and requires further extension [5]. Along with the elliptic functions 
and meromorphic elliptic integrals, the logarithmic integral should be involved to the complete 
theory [5] . Yet another aspect needs to be kept in mind is that we get in (4) the even dimensional 
vector field and may compare it locally with a straightened one having (always and at least) 
the canonical symplectic structure [14]. Actually, this program will be realized in this and next 
sections. 

The functional form of solution to (4) follows from the widely known elliptic solutions of 
the Euler top followed by supplementing the meromorphic Legendre integral E (see [1, 18] and 
appendix for notation and details); 


X =—kasn{at + £-,k), y = ia dn(at-|- e; /c), z = \.kacii{at + £;k), 

^ = —k‘^a^ J sn^(Q!S -|- e; k)ds = K aZ{at + £-,k) — a^t, (5) 

where {k, a, K-,£) are free constants. However, we are about to get some invariant constructions 
(Poisson bracket etc) and therefore one should transform this solution into the integrals of motion 
form. 

Let (I, J, K) denote the three independent functions of (x, y, z, ^) being constant on solutions 
of (4) and let AA be a function of the same variable set evolving linearly in time. Then the 
locally equivalent form to the system under consideration is 


j\/ = l, 7 = 0, j = 0, k = o. 


( 6 ) 


We shall informally name the {J\f; I, J, K)-set as ‘straightened’ variables (tt^, ..., tt"^). Clearly, 
the pair (/, J) is an equivalent of (a, k) and M is defined up to a function of (I, J, K), that is 
the constant e. To put the scheme differently we are seeking for a formal point transformation 
(x, y, z, ^) (vr^,..., TT^) between the original and straightened variable set in both directions. 

Then some routine computations with elliptic integrals yield the sought-for change: 


X = —I sn ( JAf] 


y = i J dn ( JAA; 


C = K + JZ( - J^M, 


M = 


kx^ — z^ 


^x2 — y2 y ' \/x^ — y^ j 

^ e + \/a :2 - y 2 J F 


yjx^ — z^ 


\/x2 — ^2 ’ -^x2 — y2 


z = \I cn 



k = x‘^-z‘^, J^ = x^-y^, 

^ ^ X \/x2 — 2;2 ^ 

y-\/x2 — 2;2 ’ J J 


Hence the system (4) has only two independent single-valued (rational) integrals I and J and 
the Hamilton function H may be chosen only as their combination: H = H(I, J). Differential 
computations with all the functions (sn, cn, dn, Z, F, E) have been completely developed; see 
the well-known Jacobi calculus [1, 18] accompanied with formulas (12)-(14). 
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4 Poisson bracket(s) 

The Hamilton function and Poisson bracket H appear always in a pair and the former is subjected 
to our free assignment. On the other hand, every vector field is locally equivalent to the simplest 
one, like (6), and its Poisson structure may be assigned among the constant symplectic ones. 
For comparison with the standard Lie algebraic case (3) let us consider the linear in /^, case 
aP + bP: 

H = a(x^ — -|- b[x‘^ — y^) 

with some constants a, b. 


The method 

The scheme of finding the bracket is based on the fact that under fixed dimension any bracket 
for any hamiltonian system is locally equivalent to any other bracket. Hence we may take as 
an ansatz some ‘known and good, straightened’ bracket K-^^(7r) for known ‘simple’ vector field 
p _ and seek for its transformation to the initial field coordinates (x^). The x-coordina- 

tes H^^(x) of the bracket tensor is then recalculated by the standard transformation law 

It should be noted here that ansatz for K(7r) may contain some freedom which is not restricted by 
the constant parameters (canonical constant bracket K) because the only condition for a bracket, 
apart from antisymmetry, is the Jacobi identity. In terms of inverse matrix-tensor K = the 
condition turns out to be a Maxwell (linear) equations first half dK = 0 which, clearly, has 
infinitely many solutions. 

To continue further quantization we have to produce at least one nontrivial solution. Some 
experimental computations show that there exists an unusual bracket 
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0 
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n 




( 8 ) 


and its determinant is not a zero in the generic case: det H = (x^ — Interestingly enough, 
the upper-left (3x3)-submatrix of H (boxed) is also a bracket and provides yet another, rather 
nonstandard, bracket for Euler’s top; it is related to the e(l, l)-algebra. For reference we display 
also the commuting vector field corresponding to the integral P = — x^: 

x = 0, y = 0, i = 0, i = 2{y‘^-x^), {H,/} = 0. 


Remark 1. Informally, we may conclude that knowing the bracket tensor is an equivalent 
to knowing the differential computation rules for the functions realizing complete system of 
integrals for the dynamics. Appearance/introduction of a specific form to the functions may be 
thought of as a reduction of the full system dimension to a smaller number of functions; the 
absence of integrals means thus just a generic - nonintegrable or not yet integrated - case. If we 

have differential function calculus in one direction, say, (x) --■> (vr), then calculus in a reverse 
direction (vr) (x) is automatically induced. 
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Example 1. Consider the planar vector fields 
x = A{x,y), y = B{x,y), 

where A and B are some polynomials in {x,y) over C. Recently there was developed and 
implemented an effective algorithm [3] for computing the rational integral H = R(x, y) to these 
equations. In this connection P. Olver put a question (October 2014): what about Hamiltonian 
structure for these systems? As we have seen from scheme (7), this question has a constructive 
answer. 

First, let us straighten the field under question. One of the = (A?, i7)-variables 

we’ve already had; this is the rational 77-Hamiltonian R{x, y) coming from the algorithm men¬ 
tioned above. The second variable - momentum conjugated to 77 - is derived from equations 
determining trajectory of the system. Indeed, 

R{x, y) = const 

and, hence, variable N is an integral 

M=[ =: )X{x\R{x,y)) 

J A{z,w) 

and this integral is an Abelian one with respect to algebraic irrationality 

72 ( 2 ;, tc) = const. (9) 

Here const) is a function representation to the integral as function of the x-variable and 
a constant being a modulus of the algebraic curve (9). Differential computations for this Abelian 
integral should be developed independently; including, which is an important and nontrivial 
point (see Remark 1), differentiation with respect to ‘const’ of the curve (9). Since H is an 
Abelian integral its differential calculus will possibly involve Abelian integrals of other kinds 
as well as rational functions of the pair (x,?/) [2]. Anyway, development of the calculus is 
a technical and solvable problem because derivative of the Abelian integral in a parameter is 
again an integral belonging to the same algebraic curve. 

Now, we assign the function 72(x, y) as a Hamiltonian 77 = R{x, y) and take the simplest 
canonical local (2x2)-bracket K = (_? g) straightened vector field 

^ = 1, 77 = 0. 

Then the sought-for (x, y)-coordinates Q{x,y) of the bracket is calculated explicitly by (7). 
Suffice it to display only the D^^-component: 

^ = nVR, = ^yR ■ - V^72 • (10) 

Rationality or the global single-valuedness of the found bracket tensor is established (if any) on 
the base of differentiation rules for the Abelian integrals that arise. Then it remains only to 
reassign the Hamiltonian 77 —)• /(77) for a simplification of the final answer. 

Let us consider a further specific subexample: the vector field 

X = y, y = 6x^ 



7 dx 

^ J A{x,y)' 


defining the equianharmonic elliptic Weierstrass function x = p(t; 0, o) and its derivative y = 
p'(t;0, a) [18]. According to the scheme above, we have 


77 = l(y2_4x3)^ 



V 4s^ — const 


=> A? = 2l(x; 0,4x^ — 
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where 2l(x; 52) ffs) ^ notation for a holomorphic integral in the Weierstrass elliptic theory [18]. 
Although differential calculus with Weierstrassian integrals is also absent in standard literature, 
it can be developed in the same manner as Legendrian formulas (12), (13). Computation (10) 
with rational function R{x,y) and transcendental one 2l(x;0,— y^) yields: = 1. This 

(global) unity corresponds exactly to the standard canonical bracket Q = (_? g). It should be 
emphasized here that the seed ‘straightened’ bracket K is not well-defined but is just a ‘local 
unity’ since J\f = 21(x; 0, — y^) is a multivalued function. Meanwhile the (x, y)-coordinates of 

the bracket are the global well-defined ‘object unity’. Clearly, the (2x2)-case under consideration 
has an additional convenience. Reassignment of the Hamilton function exhausts all the possible 
brackets and there is no need to investigate nonconstant ‘straightened’ brackets K{N,H). 

Example 2. In a slightly different but equivalent manner the Hamilton r-dynamics was found 
in work [6] for the elliptic 0(O|r)-constants; the very ^-coefficients of system (1). 

Example 3 (exercise). Extend the method into odd dimensional cases and reproduce the stan¬ 
dard Poisson SO(3)-description (3). 

5 Quantization: operators and spectrum 

In this section we propose a complete canonical quantization procedure for the model (8) de¬ 
scribed above. As a byproduct, we get also some nonstandard quantization for the classical 
Euler equations (3). 

As usual, we should, according to the correspondence principle 
{x,y} --4 [x,y], 

find representations for the phase observable operators {x,y,z,^) --■> (x,y,£,^), where bracket 
{x,y} is given by the H-matrix (8). Taking into account that Euler’s (3x3)-subsystem (3) is 
realized by the 1st order differential operators with linear coefficients (the classical momentum 
operators [10]), we can try the same structure ansatz for our (4x4)-extension 

(phase variable operator) = (linear function) x ^variable + (linear function). 

Note, the standard generators dx, dy, xdy + ydx of the aforementioned e(l, l)-algebra are also 
nicely fit this ansatz. Such a solution - combination of dilatations, rotations, and (9-operators - 
does exist and we find 

X := X, y := y, T := -xdy - ydx, ? := xdx + ydy (11) 

(first pointed out to the author by P. Kazinsky). Commuting relations 

[x, z] = y, [y, £] = X, [?, x] = x, [?, y] = y 

are the same as the bracket ones for base variables (x,y,z,^), i.e., (8): 

{x,z}=y, {y,z} = x, {^,x} = x, {?,y} = y 

(all the absent pairs are zeroes). For simplicity we omit here the physical multipliers like ih in 
front of commutators. From (11) it follows immediately that quantum Hamiltonian is repre¬ 
sented by the operator 

:= ^(z^ - x^) = ^{(xdy + ydxf - x^} 
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and we readily check that quantum equations of motion (Heisenberg’s dynamics) also holds with 
respect to the canonical Weyl ordering 

X = ^(pz + ^) = [x,Jif], 

y = 

^ = ^i^ + yx) = [z,J^], 

i = -^(xx + xx) = 

Availability of a (complex) rotation operator xdy + ydx in Jif prompts us to make a variable 
change to the radius-angle (r, 7 )-variables 

x = 2 rcos(| 7 ), y = 2 irsin(| 7 ), 

whereupon we get z = —2d-y and the following equivalent of the eigen-value problem 

JF'if = =-\i^E+ r^)^. 

Variable r (commuting observable I) has been separated, i.e., became a parameter, and we 
may finally normalize the spectral problem into the Schrodinger equation with a trigonometric 
cos-type potential 

= (Acos 7 - E)^ 

and two extra parameters E and A. This is the very well-known equation - Mathieu equa¬ 
tion [1, 18] - and excellent book on its theory has hitherto remained the MacLachlan treatise [12]. 

Until now, we were dealing with integrable but complex equations. Their general integrability 
is related essentially to the complex phase space but for physical motivation we need a real 
reduction of the space as long as the reduction was a non-contradictory one; see also some 
arguments on reality and complexification at the end of § 1 in [15]. A rational question at this 
point is: Whether this model is realistic? In other words, is it possible to ‘adjust’ our formal 
(yet complexified) quantities in system (4) to be real physical observables and their operators 
to be (real) self-adjoining ones? The question has a positive answer. Indeed, we want to have 
real H and, hence, x^,z^, 7 , G M under 7 ^ 0 (detU 7 ^ 0). These requirements can be 

satisfied if we put x, ^ G M, y G iM, 2 ; G iM and 7 , r G M. So, changing 2 —iz, y —)• iy we keep 
the system (4) to be pure real and Hamiltonian 

--4 + x^) 

then becomes not only the real operator but also a lower bounded one. It immediately follows 
that background state and quantum mechanics do exist and therefore this Hamiltonian may be 
‘not a bad’ candidate for the realistic one with a physical spectrum. 

Mention should be made of relations between proposed scheme and some closely related 
models. On the one hand, the Mathieu equation can be treated as a quantum physical 1-di- 
mensional pendulum and its limit A cos 7 = A(l — ^ 7 ^-!-•••) reproduces locally a quadratic 
oscillator with an infinite discrete spectrum. Parameter A may play here a role of the potential 
well depth. The model (4) does thus inherit properties of these limit cases (a dimension reduction 
4--->2 and approximation cos 7 ---> 7 ^) and, on the other hand, is an alternative to the standard 
SO(3)-quantization of a rotator [10]. It is interesting to observe that quantum statements 
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E = const 


A = const 


0 14 9 16 E 


Figure 1. The gap structure spectrum for the Mathieu equation 4'" = (Acos-f — E)'^/. Shaded domains 
are the spectral lacunae (‘forbidden’ zones). 


in both this quantization scheme and proposed new one are determined by the complete set 
of compatible/commuting observables which are classically evolving as meromorphic functions 
of t; not as meromorphic/logarithmic elliptic integrals. Recall, meromorphic (elliptic) functions 
are the single-valued (well-defined) functions on torus, whereas integrals, in general, are not. 
Indeed, the standard momentum square and for the rigid Euler’s body are combinations 
of sn-, cn-, and dn-functions [10]. In the model under consideration the variable which is 
a meromorphic elliptic integral (5), does not get to the commuting pair of observables and 
all the other phase variables are also the meromorphic functions on a torus. We observe in 
passing that integral structure of the classical solutions to this new model is much simpler 
than the standard one. One of the Euler classical position (-0, 9, (^)-angles of a body - the 
angle (p{t) - is a cumbersome combination of meromorphic and logarithmic elliptic integrals; 
see [9, Chapter VI, § 37, equations (37.17)”(37.20)] and [11, p. 139]. Meanwhile the only variable 
of an integral type in the model (4) - variable ^ - is just a simple Z-function or, which is the 
same, some 0i/0-combination. One may be noted in this connection that complete quantization 
procedure of theta functions will perhaps inevitably involve, the holomorphic integrals excepted, 
the full set of Abelian objects on Abelian variety: meromorphic functions, meromorphic and 
logarithmic integrals, and, which is the most nontrivial point, the 6 itself as integral of an 
Abelian meromorphic integral. Such a nature of the theta-function is explained in [5] . 

Einally, Eig. 1 exhibits the typical spectra structure for the Mathieu equation [13] in both 
the quantum numbers E and A. The number of lacunae (colored in red) is infinite and bold 
green dots on the horizontal E-axis is the sequence 0^, 1^, 2^, .... It is a gap boundary set 
of the limit case A —)• 0 following from the fact that the zero/constant potential has pure 27r- 
periodic solutions when spectral parameter is an integer perfect square E = . Exhaustive 

theory, numerous applications, and wealth of numeric information to the Eig. 1 can be found 
in [8, 12, 19]. Note that number of gaps in A-direction is also infinite as in the E-axis; for proof 
see [19, Chapter 7, § 3.5]. 
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A Appendix. Differential calculus of Legendre integrals 


The standard Legendrian elliptic integrals [1] are the holomorphic and meromorphic integrals 

ds {\ — k‘^s^^ds 


F{x-k) := / 
JO 


0 -^(1 - s‘^){l - k'^s^y 


E{x;k) := [ 
Jo 


0 v^(l - s^)(l - k'^s'^) 


(1st and 2nd kind respectively) and the logarithmic integral (3rd kind) 



1 ds 

1 - ^(1 _ s2)(l - fc2s2) • 


The latter has a parameter a. 

The triple (-F, E, IIq) forms a linear vector space and the theory above requires complete set 
of its partial differentiations dx, dk, and Oq. Somewhat surprising facet is the fact that no self- 
contained differential calculus of these objects with respect to all the arguments {x, k, a) seems 
to have hitherto been tabulated in handbooks; see, e.g., [1, 7, 16]. Bearing in mind a remark 
about logarithmic integral at the beginning of Section 3, we make up this deficiency. 

Denote := (1 — x^)(l — Then one can show that the three objects {E, E, n„) have 

the following closed rules of differential computations 


dE 

dx 
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dE 

dk 


F 
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xy 


dE 0^1 dE F E 


dx " 'y 

dYia _ 1 1 

dx 1 — y ’ 
-k 


k{k?‘ — 1) — 1 1 — /c^x^ ’ 

' _ F E 

dk “ ^ ¥ 


k'^xy 


dk (/c^ — q/)(/c2 _ 
1 


da 


2a(a — l)(k^ — a) 


E + (k^ - l)na - 
(k^ — a)F + aE — (k^ — — - 


a^xy 
— ax^ 


( 12 ) 


(13) 


Note that these equations are compatible; the cross derivatives like 
d dE _ d dE d dUa _ d 

dx dk dk dx ’ da dk dk da ’ 

turn out to be identities. If we introduce also a Jacobi analog of Weierstrassian ("-function, i.e., 
u-representation 


Z(u; k) := £l(sn(u; fc); k) 


of the ii-integral, then it is differentiated as follows 

^ = 1 —/c^sn^, ^ —(Zcn^ — (fc^ — l) usn^ — sncndn). (14) 

au dk k^ — 1 

If required, partial /c-derivatives of functions sn, cn, dn(u; k) are calculated by formulae which 
are known [11, pp. 82-83]; we do not display them here. 
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